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We present a method for parameterizing the charged particle multiphcity and charged hadron 
spectra from heavy ion data with a simple, Glauber-inspired, model. The basis of this model is 
derived from the observation of leading hadrons in pp collisions and the number of interactions 
calculated by a Glauber model. Singly hit and multiply hit nucleons are treated as different sub- 
components of the same collision. With this scheme, for fixed sub-component yields, we find that the 
multiplicity and charged hadron spectra can be reproduced, without the need for large suppression 
at high-pT. Suppression is still observed, in the low- to intermediate-pr region, but this is confined 
to suppression of surface emission partons. At high-pr, the suppression disappears entirely, leaving 
only a possible Cronin-type enhancement. The suppression observed in 200 GeV collisions is found 
to be the same for 62.4 GeV data. 
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INTRODUCTION 



The comparison of hadronic collider {pp collisions) and electron collider (e^ + e" collisions) data have shown a clear 
discrepancy in the produced particle multiplicity for collisions measured at the same center-of-mass energy. This has 
long been explained by the observation that hadronic interactions can be decomposed in terms of "collision products" 
and "leading, hadrons" - with the latter carrying away (on average) half of the total energy available for particle 
production [1|. It is found that one can unify these data by shifting the center-of-mass energy of the pp collisions to 
an effective center-of-mass energy, which is half of the original collision energy, empirically determined. Superimposing 
central heavy ion data onto this agrees more closely with the e"*" -I- e~ data than the hadronic collision (pp) data 0. 
Combining all these observations, one can surmise that collisions in which all "constituents" of the proton are used 
(i.e. no leading hadron) will produce higher (e+ + e^-like) multiplicities. 

The model presented here, originally described in Ref. [3], uses these ideas as a premise to divide the participant 
region in heavy ion collisions into two distinct regions; these are noted as "mono" and "multi" . A Glauber model [3| 
is used to calculate the number of participants and the number of times each nucleon is struck by another from the 
opposing nucleus. Mono refers to nucleons which undergo a single hit - similar to pp interactions, which presumably 
can "liberate" a leading hadron. Multi refers to nucleons which are multiply struck (by two or more) nucleons, 
however, these are only counted once and thus do not accumulate as quickly as the number of binary collisions. 
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FIG. 1. (color online) Predicted relative positions of mono- (upper) and multi-collisions (lower). The distributions are 
based on a Glauber calculation at impact parameter 6=2, 5, 8,12 fm (A'part ~ 359, 261, 140, 24, A'^.^ono ~ 28, 38, 33, 12, and 
A^muiti ~ 331, 223, 107, 12). In this picture, the reaction plane is always along the s-axis. 



This phenomenological approach is similar to the successful core/corona approach of EPoS [5|- The main technical 
difference surrounds the expected position with respect to the reaction plane, of the mono (or corona) nucleons. The 
current model expects that the peak of mono nucleons to be along the reaction plane (see Fig. [1} , whereas in the corona 
from EPoS is expected to peak 90" away [H]- A more recent 'toy model' implementation by the EPoS authors @, 
however, does use the same "Glauber-method" to determine core (multi) and corona (mono) regions to describe the 
strangeness enhancement observed in Au-|-Au and Cu-|-Cu data. In that implementation no interactions between the 
core and corona regions is assumed, each "corona" contributes one-half of a minimum bias interaction. The authors 
note that fixing the mono (corona) is not necessary. However, although they undergo different assumptions, both can 
provide a good description of the data. Where these models differ is in the extent to which they can be (or at least 
have been) applied and the possible origin of the two components, or regions, of the collision. By using the model 
presented here, it is possible to not only describe with success \ow-pT phenomenon observed in heavy ion collisions 
but also high-pT phenomena. This can lead to hints as to the origin of disappearing back-to-back jets and their 
reappearance at low-pT- 
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FIG. 2. The relative number of mono and multi collisions versus the number of participants. Black (grey) circles repre- 
sent Au+Au (Cu+Cu) collisions at 200 GeV (ctjvjv = 42 mb); closed and open symbols represent mono and multi collisions 
respectively. The values for 19.6 GeV (crjvjv = 33mb) Au+Au data are depicted as lines for comparison. 

THE MODEL 

The implementation of this model is very simple. A Glauber model is used to count the number of times each 
nucleon is struck (the Glauber model implementation is similar to Ref. 0). For each Monte-Carlo "collision" the 
nucleons are sorted into three categories: not hit (spectators), singly hit nucleons (A^mono), and multiply hit nucleons 
(TVmuiti)- Spectators are not considered in this model. After many trials, and at all possible impact parameters, 
the average number of mono- and multi- assigned nucleons are calculated. Figure [2] illustrates the mean A^mono and 
A^muiti found for Au+Au (black circles) and Cu+Cu (grey) collisions at ^snn = 200 GeV. N^ono and A^muiti versus 
centrality are not found to appreciably change with collision energy, as shown by lines depicting A^mono and A^muiti for 
19.6 GeV collisions. 

After calculating the number of mono and multi participants for a given centrality class it is possible to form the 
multiplicity (and spectra, etc.) corresponding to that data using the formula: 

Ya u+Au — -^mono^^iono ^~ -^multi^nulti (1) 

Yau+Au represents the modeled Au+Au data, ymono is the yield expected from singly interacting nucleons (ostensibly 
minimum bias pp data) and Ymuiti is the yield expected from the multi underlying distribution. Fmono and Fmuiti are 
fixed for a given data sample and do not change versus centrality. Y itself represents any distribution, for example 
multiplicity, charged hadron spectra, kinetic freeze-out temperature, or any other observable. 

In this paper, we first explore the possible origin of the underlying multi distribution which could be from one of 
many different sources such as non-single diffractive (NSD) events (unlikely as these may still have leading hadrons), 
high-multiplicity events (for example greater than (N^^) /2.), high-p-r data (pt>~1-8 to represent the close-packedness 
of the multi nucleons) , a gluon-only source distribution. We use each of these different hypotheses to test the systematic 
dependence - although not to choose a preferential particle production mechanism but rather to provide a proof of 
principle that one can reasonably reproduce the particle distributions in heavy ion collisions from these simple sources. 
We note that a single distribution for mono (pp minimum bias) and a single distribution for multi is assumed, which is 
sufficient to reproduce all the data versus centrality. The centrality dependence observed in the data is then entirely 
due to the interplay between iVmono and iVmuiti, Fig. [21 and not from a changing or modified underlying distribution. 

Secondly, as a more rigorous test, fits to the data to extract the underlying mono and multi distributions will be 
made. In fact this is necessary to analyze more exotic data than multiplicity and charged hadron spectra. In this way 
we decouple ourselves from (for example) Pythia Q to describe the detailed identified particle spectra. By using 
Eqn. [T] and the measured data to simultaneously extract the underlying distributions, both the modification to the 
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pp minimum bias collisions at the surface (mono) as well as the true underlying multi distribution can be studied. 
As Glauber has no recourse to produce particle distributions, this model has a purely phenomenological approach to 
describe the heavy ion experimental data. 

One important distinction to note for the current model versus many other two-component models is the absence of 
"collision scaling" . Other two-component models assert that one component factorizes with the number of participants 
(usually associated with soft particle production), and the second with the number of nucleon-nucleon interactions 
(-^coU - representing hard collisions) within the collision [§- 12|. In the current discussion, we assert that the number 



of collisions is not the correct scaling variable. If, as is hypothesized, upon collision the system becomes a hot, 
dense QCD matter then it may become difficult to reliably count individual nucleon-nucleon collisions. The currently 
accepted approach used to validate the number of collisions is to measure the spectrum of direct-photons in A+A 
collisions, relative to that in pp interactions. As the direct photons are not expected to interact with the medium, it is 



assumed that the measured Raa = 1 validates the use of A^coU uM ' illustrating that the number of originally produced 
photons are consistent with pp. In more current data, however, the previously observed iVcou scaling at pt ~ 6 GeV/c 
is not seen at higher momenta IJ] , even though punch-through high-p^ hadrons are seen to escape the medium [15| . 



Here, we set this A'coU scaling assumption aside. The number of mono and multi collisions in the current model is 
less susceptible to the assumptions needed in order to count the number of collisions, where our number reflects an 
interaction volume (multi) and simply the number of expected nucleon-nucleon scatterings on the surface (mono). 

For a low number of participants (A'part<25), the number of mono collisions exceeds that of multi, and vice- versa for 
a high A^part- However, as shown in Fig. [2l the number of mono does not fall to zero even for the most central events 
where approximately 7% of the participants undergo only a single interaction. An important consequence in these 
calculations is that the sum of mono and multi, by definition, is fixed to be A^part- This restricts the growth of multi, 
unlike the rapidly growing number of binary collisions, and results in multi (and mono) being more closely associated 
to a "participant" variable. In calculating the mono and multi, it is found that there is a strong dependence on the 
collision species, in contrast to the number of binary collisions (A"cou) which is essentially the same at the same A^part- 
Similarly, only a weak dependence on the collision energy (inelastic pp cross-section) is observed over the energy range 
of the RHIC data, whereas a strong binary collision dependence on A"coii is noted. 

MULTIPLICITY DISTRIBUTIONS 

As the premise of this model uses the charged particle multiplicity to argue for distinct mono/multi regions of 
the collision volume it is natural to begin with a discussion of that experimental data. Here, one can simply use 
the measured pp distribution as the mono and a second underlying distribution (for example taken from Pythia) to 
represent the multi. We combine the mono/multi distributions (and similarly for the next section) as: 

.A^Au-l-Au JTYmono JTymulti 
; — jVniono ; H JVniulti ; i^j 

0,7] dr] drj 

where Amono and Amuiti distributions for Au+Au collisions are given in Fig. [51 One can hypothesize that dN^°''^° / drj 
be dN^^/dr] (i.e. precisely from minimum bias pp interactions) or it can be a free parameter in the fit. For simplicity 
in this first case, we assume the former; later, the latter is used to test for any modification of the underlying mono 
distributions. 

Figure |3] shows the result of using this formalism with several hypothesized multi multiplicity distributions for 
Au-I-Au collisions at ^snn = 200 GeV. For the shown centrality bins from PHOBOS |16jl, one can see that the model 
provides a reasonable description of the data. Remember that the same underlying distributions are used for both 
panels, only A^mono and A^muiti are changed. This is not too surprising as the basic premise of this model is to 
reproduce the centrality dependence of the charged particle multiplicity in the context of an absence of leading 
hadrons in the multi region. From this, judging the model's success is really only dependent on the choice of the 
underlying multi distribution. Three simple pp multiplicity-like distributions are used for the multi component, for 
comparison. The pp collisions at ^/sm^ — 200 GeV minimum bias distribution is used (closed circles) which is similar 



to the "Wounded Nucleon Model" [17| approach and yields a multiplicity which is ~30% lower than the measured 
Au-|-Au data. Rejecting the single diffractive events yields a higher multiplicity (open circles) but is still significantly 
lower than the Au-|-Au distributions. A more radical approach is to only choose the highest multiplicity pp events, 
in this case more than the average (>24) (squares), which does achieve the desired multiplicity. One can understand 
this approach in terms of the origins of the model - as higher multiplicity events typically have lower energy leading 
hadrons, thus releasing more energy to particle production. 
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FIG. 3. (color online) Comparison between measured multiplicity versus rj (PHOBOS ^snn ~ 200 GeV Au+Au data 16;]) 
(band) and model predictions with various candidates for the underlying "multi" dN^^^^/drj distribution as noted in the 
legend. 0-6% central (left) and mid-central 35-45% (right) data are shown. In each model representation, the "mono" is fixed 
to the inelastic pp distribution at the same energy, i.e. dN'^'^° / dr) = dN^^ / d-q. Note the t/-axis scale difference in each figure. 



Changing the reference pp multiplicity is not the only way to alter the resultant multiplicity distribution. Also 
shown are a subset of events from the minimum bias pp cross-section which exclusively selects out gluon interactions 
{99 ~^ 99 or 99 ~^ 99) with a large momentum transfer {pT > 1.8 ± O.lGeV/c) (lines). (Note that the same 
multiplicity enhancement is also found by using all (quark and gluon) interactions with px > 1.8GeV/c.) Such an 
approach may represent the close-packedness of the nucleons which may facilitate a large momentum transfer. The 
different minimum pT values used as the multi component (denoted as dotted and dashed lines in Fig. [3]) show the 
sensitivity of the choice of px in the final distribution (dN^f^^^^ /dr] of Eqn. [2); for multiplicity this is small. 

Rather than assuming a given Fmono and/or Fmuiti multiplicity distribution, Eqn. [T] can be used to extract these 
two underlying distributions from the data. Two bins in centrality (preferably far apart in centrality) are used to 
simultaneously extract Fmono and Fmuiti- This is referred to as the pair fit method. After finding the underlying 
distributions, these are fixed and all centrality bins are recombined using the A^mono and A'muiti as in Eqn. [T] To 
test the sensitivity to the bins chosen, several combinations are made. The top panels of Fig. |3] show the PHOBOS 
data (grey bands) with the fit distributions (solid lines). Panels (a,b), (c,d), (e,f), (g,h) in Fig. |31 show data for 
Au-hAu coUisions at 200 [11], 130 [Til, 62.4 [ll], and 19.6 GeV [20| respectively. In each case, the model uses the 
most central 0-6% and 35-45% (mid-central) data to extract the underlying mono and multi distributions. Using 
other bin combinations does not significantly alter the resultant distributions. The lower panel of Fig. U] shows the 
underlying Yinono (dashed line) and ymuiti (solid line) distributions from the fit. The grey and open bands represent 
the systematic error of the data. For comparison, the gluon-only distribution selected from Pythia (dot-dashed lines 
for the pt > 1.8 GeV/c distribution only at 200 GeV) along with the minimum bias at the same center-of-mass energy 
(circles) are also shown. The underlying multi distribution is considered "stable" in the sense that using different 
centrality pairs does not change appreciably the extracted Fmuiti distribution. Mono, however, is less stable when 
using the pair fit method and is found to be more dependent on the choice of centrality bin. Several interesting 
artifacts of the fit are observed. First, the extracted Ymuiti distribution is close to that which would be obtained by 
using a gluon-only pp distribution in Eqn. [5J but with more yield at mid-rapidity and also narrower. Second, the 
mono distribution is considerably lower at mid-rapidity than the measured minimum bias pp collisions pll . [22j , as 
shown by the comparison of the mono-collisions from the pair fit to minimum bias pp in Fig. |4l panels (b), (f), and 
(h). This is possibly hinting that a modification exists, due to an opaque multi region. A more detailed description 
of this modification is discussed in the next section. The final notable observation is that the yield of the extracted 
mono interactions is higher at forward rapidities, as compared to the minimum bias pp. This could possibly point to 
an influence of spectators in the measured data. To reiterate, Vmuiti represents the multiplicity per participant pair of 
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the multi region, which is consistent with the expected yield from a gluon-dominated collisional center. The core is 
found to be opaque, with some of the multiphcity from the surrounding mono colhsions absorbed. Importantly, even 
though the size of the central multi region changes, the yield per participant does not change. 



(a) Au+Au 200 GeV (dala+model) 




(b) Au+Au 200 GeV (model constituents) 
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(c) Au+Au 130 GeV (data+model) 
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Multi-Collisions 




(e) Au+Au 62.4 GeV (data+model) 
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FIG. 4. (color online) PHOBOS multiplicity data (bands) from Au+Au collisions compared to the multiplicity distributions 
derived from the model, lines, using the pair fit method (upper panels). Panels (a,b), (c,d), (d,e), (f,g) show data for Au+Au 
collisions at 200 TS], 130 ^T^, 62.4 and 19.6 GeV ^20] respectively. The lower panels show the underlying Ymono (dashed 
line/open band) and Ymuiti (solid line/filled band). For reference, the minimum bias pp data [2ll. [2^ (green closed circles) 
and gluon only distribution (green dot-dashed line) are shown (the latter for 200 GeV only). The open (closed) star symbols 
represent the result from the least-^^ method to determine the underlying mono (multi) from mid-rapidity data, see Fig. \5\ 
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FIG. 5. PHOBOS yj, I2a425l. |27|1 and ALICE ^2g] multiplicity data at mid-rapidity compared to that derived from the fit, 
using the least-x^ method, for Au+Au (Pb+Pb) data in the left panel and Cu+Cu data in the right panel. Fits are shown for 
RHIC energies ysiS" = 200GeV (lightest grey band), 130 GeV (light grey Au+Au only)), 62.4 GeV (dark grey), and 19.6 GeV 
(22.4 GeV for Cu+Cu) (darkest grey). The ALICE data (2.76 TeV) are shown as a no-fill outline. For comparison, in the right 
panel, the expected Cu+Cu data are shown using the ymono and Vmuiti yields found from the fits to the Au+Au collision data. 

There is another way to extract the ymono and Ymuiti yields from the data. Figure [5] (left panel) shows the mid- 
rapidity Au+Au data from PHOBOS fll, [ij-list and Pb+Pb data from ALICE 126*1 data (bands) together with a 
least-x^ fit (lines) using the A'^mono and A^muiti distributions from Fig. [21 The right panel of Fig. [S] shows the Cu+Cu 
data from PHOBOS [27| along with fits using that data and the extracted parameters from the fit to the Au+Au 
collision data at the same energy. (Note that the 22.4 GeV Au+Au reference is scaled by 1.04 to account for the 
small difference in collision energy.) This least-x^ fit is a different approach than the pair fit method used above. 
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FIG. 6. (color online) Mean mono and multi radii - estimated from the Glauber model. The left (right) panel shows the radii 
along the (90° to) reaction plane (the line along the nuclei centers - x- (j/-) axis in Fig. [ij . The solid line shows the multi 
radius (defined as 10% of the peak). The different dashed lines represent the mono peak position (dash) and the outer/inner 
radius, defined as 10% of the peak yield. Negative radii (in the inner radius) represents the point at which the mono are no 
longer distinctly separated by the multi - i.e. mono collisions occur at every point throughout the collision area. 



TABLE I. Mono and multi yields (and relative to data pp inelastic yields [2ll. [2^. [2^. [25j) as derived from the least-x^ fit to 
the mid-rapidity data from PHOBOS (Au+Au and Cu+Cu) and ALICE (Pb+Pb). 



Energy 
(GeV) 


dNl 
yield 


Au- 
yield/pp 


hAu 
yield 


r'^'/dn 

yield /pp 


yield 


Cu- 
yield /pp 


hCu 
yield 


r'''/dn 

yield/pp 


2760 
200 
130 
62.4 
19.6 


1.35 
1.49 
1.04 
0.86 
0.99 


0.65 
0.46 
0.45 
0.78 


8.55 
4.01 
3.56 
2.88 
1.94 


1.75 
1.74 
1.50 
1.52 


2.07 

1.58 
1.53 


0.91 

0.82 
1.20 


3.83 

2.75 
1.98 


1.67 

1.43 
1.56 



This enables the full centrality dependence of the mid-rapidity results to be simultaneously utilized and reduces the 
dependence of the results from the single pair choice. The resultant least-x^ fit at each energy well represents the 
data, even over two orders of magnitude of collision energy. The derived mid-rapidity mono and multi multiplicities 
are summarized in Table HI For comparison, these yields are shown as a star symbol to compared to the pair fit 
method in Fig. 21 

In an attempt to reconcile the difference between the minimum bias pp data and the derived underlying mono- 
distribution, we compare the detailed spacial positions of the mono and multi interactions, similar to those shown 
in Fig. [TJ Figure [6] shows the mean radius of the multi (solid line) and mono (dot and dot-dashed lines) defined at 
10% of the peak. The center of the mono (peak) is shown as the dashed line, which coincides with the multi radius. 
The left (right) figures show the radii along (90° to) the reaction plane. Clearly, half of the mono distribution resides 
within the multi distribution. One could well imagine that if the multi represents an opaque medium, then much of 
the mono inside the multi region could be absorbed (or at least modified/suppressed). Similarly, those outside the 
multi region may be suppressed should the particles be headed toward the multi region, an eclipse- type effect. 

One can test the model by using the Cu-|-Cu data to extract new l^ono and Fmuiti yields from the data and 
compare those to the Au+Au underlying mono and multi distributions. Figure [7] shows a similar analysis as described 
above for the pair fit method, this time with Cu-f Cu collisions at -^/snn = 200, 62.4, and 22.4 GeV. For completeness, 
the underlying ymono and ymuiti distributions from Au+Au collisions (blue dashed lines) are also used to form the 
multiplicity data. In comparing the underlying mono and multi distributions from Cu+Cu data to those from Au+Au 
data, we find that the multi distributions are very similar and that the mono distribution in Cu+Cu more closely 
resembles that from minimum bias pp interactions. Similarly in the right panel of FigJSl the same least-x^ method 
used for Au+Au collisions is applied to the mid-rapidity Cu+Cu data from PHOBOS [27[ where the same trends are 
observed in terms of the underlying distributions. 



UNIDENTIFIED CHARGED HADRON SPECTRA 



Following the success of the analysis of the charged particle multiplicity, we apply this approach on more detailed 
data to check whether a reasonable agreement could be found. First, the unidentified charged hadron spectra versus 
transverse momentum {pt) is used. This covers a different kinematic region than the low-p^ dominated multiplicity. 
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FIG. 7. (color online) PHOBOS multiplicity data from Cu+Cu collisions [20| ] compared to the multiplicity distributions derived 
from the model, using the pair fit method (upper panels, red) and from using the underlying Ymono and yinuiti distributions from 
Au+Au collisions (blue, see Fig. [4]). Panels (a,b), (c,d), and (d,e) show data for Cu+Cu collisions at 200, 62.4, and 22.4 GeV 
respectively. The lower panels show the underlying mono (dashed line/open band) and multi (solid line/filled band) extracted 
from Cu+Cu data. The minimum bias pp data [2ll . [2^ (green circles) and the underlying mono and multi distributions from 
the Au+Au analysis (blue dashed lines) are shown for reference. The star (square) symbols represent the least-x^ fit to the 
mid-rapidity Cu+Cu (Au+Au) data, see Fig. [S] right (left) panel. 



The suppression of high-pT particles in central events is a fascinating phenomena which has been modeled by many 
different theories, which are generally decoupled from the low-p^ region. Here, the full pt range of the data is used 
with the primary goal of a qualitative description of the data. As an initial test, the ymono and l^uiti are again taken 
from Pythia, at a collision energy of 200 GeV, with mono derived from a truly minimum bias sample and multi from 
several hypothesized distributions as described above. Again, A^mono and A^muiti are from Fig. [21 

Figure [H shows the nuclear modification factor, Raa - see Eqn. [3l for Au+Au collisions at -^/snn — 200 GeV for two 
centrality bins as measured by PHENIX [28:]. It can be clearly seen that, although the model distributions do not 
match the experimental data precisely, the qualitative features of the data are reproduced. Using the pp minimum 
bias (circles) or high multiplicity (squares) simulations for the multi component of the hadron spectra do not represent 
the data enhancement in the intermediate px region with either of these as the underlying multi distribution. The 
minimum bias result, in this case, represents a simple visual scale of the number of collisions. 



Raa = 



dydpj 



coll 



dydpT 



(3) 



nNp^rt _ dydpT / .\ 

AA nj d^NPP ^ > 

dyapT 

Applying a minimum pT cut of 1.8GeV/c on the minimum bias sample (dot-dashed line) results in a clear peak 
close to the minimum cut (as expected) but over predicts the higher-py data. Removing all quark interactions, leaving 
a gluon dominated system (solid, dotted, and dashed lines) both reproduces the intermediate-pT peak and the lower 
yield at higher-py. To emphasize again, the objective of this analysis is not to find the underlying distribution, but 
to show that a simple underlying distribution could exist which can describe the features observed in heavy ion data. 

Using the same technique to fit the data as described above for the multiplicity 'pair fit', charged hadron spectra 
can be modeled using PHENIX data [1^. The fit is performed twice, once using the 0-10% and 60-70% cross-section 
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FIG. 8. (color online) Comparison between measured nuclear modification factor, Raa, versus pr (PHENIX data (band) 
and model predictions with various candidate Vmuiti distributions (colored symbols) as noted in the legend. 0-10% central (left) 
and mid-central 40-50% (right) data are shown. In each model representation, ymono is fixed to minimum bias pp from Pythia. 
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FIG. 9. Nuclear modification factor (scaled by A'^part not A'^coii) for unidentified charged hadrons from PHENIX [2g| (filled band) 
at y'sNN ~ 200 GeV. The lines are the R^^'* reconstituted from the underlying mono- and multi-distributions using the pair 
fit method; dashed (dotted) lines use centrality bins 0-10% (20-30%) and 60-70%. In each figure, the dot-dashed line represents 
the expected A^coii scaling, assuming particle production scales as pp x A'^coii. 



bins as the data pair to obtain the Fmono and Fmuiti distributions and a second pairing of 20-30% and 60-70% to test 
the sensitivity in the choice of bin pairing, see Fig. |9l The extracted underlying mono and multi distributions are 
then used to model i?^^"", Eqn|31 in the remaining centrality bins. 

We find that all multi-dominated centrality bins are well reproduced from the underlying mono and multi distribu- 
tions using any two bins; the results are not sensitive to the choice of centrality bins. In each figure, the double-dot 
dashed line represents the '-/Vcoii' scaling. For data presented as Raa, the line would be exactly at unity and the data 
yield (at high-p^) would represent the apparent suppression. Within the mono/multi framework, we find that sup- 
pression at high-pT is not necessary to describe the data. The most peripheral bins, perhaps, are not as well described. 
An explanation for this limited agreement could be that a "collectivity" which produces the multi distribution could 
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FIG. 10. (color online) Mono and multi nuclear modification factors (scaled by A^part not A'coii) using the pair fit method on 
PHENIX, PHOBOS and STAR HO] data. The left panel shows ^siS"=200GeV mono and muhi distributions, derived 
from a fit to the PHENIX data. The dashed and dotted lines show the fits from the centrality bins described in Fig. the 
outline shows the 1-a band, derived from all possible centrality bin pairs. The grey band shows equivalent PHENIX d+Au 
data [33 from minimum bias collisions for reference. The closed (open) circles represent the results from the least-^^ method 
applied to the spectra (omitting the four lowest centrality points from the fit). The right panel shows the mono and multi 
spectra using the least-x^ fit method. The closed red (open blue) circles represent the multi (mono) collisions spectra per 
participant. The solid line represents the 200 GeV pp spectra, taken from [28^ . 




FIG. 11. (color online) Comparison of the mean mono and multi distributions derived from fits to PHOBOS (solid lines), 
PHENIX (dashed), and STAR (dotted) for Au+Au collisions at 200 GeV. 



be dissipating. This could point to a possible need to modify the model to account for an additional scaling variable 
for example mono, duo and multi-^; this is outside the scope of the current paper. 

Figure[TO] (left panel) shows the underlying mono and multi nuclear modification factors for y/sjm = 200 GeV derived 
from the pair fit method; with the corresponding spectra in the right panel. The dashed and dotted lines represent 
the centrality bin pairs used in the analysis shown in Fig. |9l The red outline shows a 1-a band around the mean 
underlying multi distribution determined from an average over all possible centrality-pair combinations, except the 
most peripheral. This appears to be well constrained, illustrating only a small variance in the underlying multi 
distribution. A comparative analysis using the least-x^ method (for the yields versus centrality for each px bin) was 
also performed. (Note that the lowest four centrality points were omitted from the fit.) The results are shown as 
circles in Fig. [10] (left panel) and are in good agreement with the distributions extracted from the pair fit method. The 
extracted underlying mono distribution, however, is less constrained, possibly due to a changing surface suppression 
of mono with centrality. Nonetheless, the extracted values of the mono distributions here are in agreement with the 
ones obtained from fits to the Au+Au charged particle multiplicity: at very low pT, a suppression is observed with 
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FIG. 12. Nuclear modification factor (scaled by A'^part not A'^coii) for unidentified charged hadrons from PHOBOS [s^ (filled 
band) at -^snn ~ 62.4 GeV. The lines are the -R^^"' reconstituted from the underlying mono- and multi-distributions using 
the pair fit method; solid (dashed) lines use centrality bins 0-6% (15-25%) and 45-50%. In each figure, the dot-dashed line 
represents the expected A'^coii scaling, assuming particle production scales as pp x A'^coii- 



respect to minimum bias pp collisions. In the intermediate pt region the mono distribution falls to zero, hinting that 
the absorption (or suppression) is maximal in this kinematic region. For some fits, it is found that the suppression is 
'negative' this would represent a case when some additional suppression is present from the edge of the multi-region; 
pointing, perhaps, toward the need for a mono, duo, and multi-'^ prescription. Figure [TT] shows a comparison of the 
resuhant distributions by using PHENIX, PHOBOS 0, and STAR [30] data. The muhi distribution in this analysis 
are found to be very similar. To put the results into context, Fig. llOl shows the multi spectra (per participant) relative 
to the minimum bias pp spectra, which is found to have the same systematic dependencies as the gluon-dominated 
system found in Fig. [8l It has been seen that, through the separation into multi and mono sub-components of the 
collision, a consistent picture can be observed. At the center of this picture is a dense (perhaps gluon dominated) 
system which is opaque to the particles produced at the periphery of the collision. A gluon-dominated picture is not 
excluded by other measurements, for example, the anomalous increase of baryons, with respect to mesons, could be 
expected from a system with a higher number of gluon-jets than quark-jets. 




FIG. 13. (color online) The left panel shows the underlying mono (triangles) and multi (circles) -R^^"'' distributions for 
PHOBOS Au-(-Au data at -^snn ~ 62.4 GeV [s^]. The solid lines are meant to guide the eye. The open (closed) symbols 
represent a pair fit using centrality bins 0-6% and 45-50% (0-6% and 45-50%). The dashed lines depict the underlying mono 
and multi distributions from Au+Au collisions at 200 GeV (PHENIX). In order to better compare the multi distributions, the 
green squares represent the 62.4 GeV data rescaled to the same pp reference (i.e. assumes a 200 GeV pp reference, not 62.4 GeV). 
The right panel has the same notation as the left, except the data is for PHOBOS Cu+Cu collision data at 200 GeV [s^ ]. 



It is interesting to note that the disappearance of the away-side jet in two-particle correlations was observed in a 
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similar pT region (for the associated particles) as the mono suppression 3l|. In this kinematic region, l<pT<3 GeV/c, 
we find that the underlying mono distribution is very small which perhaps hints that jets are dominantly from the 
mono (pp-like) interactions which occur on the surface. With this assumption, the away-side is then suppressed (or 
fully absorbed) within the (multi) medium. At higher p^, the mono and multi distributions essentially become the 
same and are similar in magnitude to the measured Cronin-enhanced d-l- Au data [32| . It should be noted that in order 
to form these spectral distributions, the A^'coii scale is explicitly removed. R^a^ distributions are formed entirely with 
our participant-like parameters, thus any large-scale suppression in the data (using Raa) is actually brought about 
by scaling the distributions by iVcoii- In other words, the suppression at high-py could be entirely an artifact of using 
the wrong scaling variable (namely iVcoii)- How could A^coii be the wrong variable? Consider a simple picture. In a 
standard Glauber model calculation, we assume that each nucleon-nucleon interaction is independent and that the 
cross-section for each interaction is the same (i.e. (Tnn does not change - even after 10 collisions). Once a "collision" 
has occurred one could surmise that the cross-section is not the same, something has been lost so the cross-section 
could be diminished. Further, if one assumes that the multi component melts into a single (not discrete) system, then 
it may be impossible, in fact improbable, to manifest multiple collisions, or at least reliably count them. Note that, 
experimentally, one cannot count A^coU, yet one can directly count TVpart, or at least derive it from the data with some 
certainty. 

It has already been observed that by dividing the -y/sNN = 200 GeV charged hadron spectra by that from 62.4 GeV, 
the ratio (for a given pt) is invariant across all centrality bins measured [27j . This may already indicate a simple 
geometry scaling which factorizes in energy and centrality. We should therefore expect a similarly successful fit using 
the mono/multi approach for 62.4 GeV data. Fig. [T2]shows the same analysis using PHOBOS charged hadron spectra 
at y^snn = 62.4 GeV [33[ , where two centrality bin pairs can be used to predict all other bins. Fig |13l shows the resultant 
mono distribution, illustrating the same level of suppression as 200 GeV, whereas 62.4 GeV has a higher muhi i?^^" ■ 
In the calculation of i?^^"'*, the pp spectrum is used in the denominator to observe any differences between A+A 
and pp collisions. To more fairly compare the two spectra, the green squares represent the same 62.4 GeV data, but 
rescaled to efi^ectively use the pp 200 GeV reference. In this way, we see that the underlying multi spectrum grows 
with energy. A similar rescaling of the underlying mono distribution is not performed as for that case, the relative 
change to pp is important, not the overall spectrum. It is interesting to note that -R^X'^' mono-collisions is the 
same for 200 GeV and 62.4 GeV. This indicates that the level of suppression does not depend on collision energy (at 
least for these RHIC energies). 
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FIG. 14. Nuclear modification factor (scaled by A^part not A'coii) for unidentified charged hadrons from PHOBOS (filled band) 
in Cu+Cu collisions at ^snn ~ 200 GeV [s^l and from the model fit using the pair fit method (lines). For the soHd Unes, the 
centrahty bins 0-6% and 45-50% were used, dashed fines used centrality bins 15-25% and 45-50%. 

In the multiplicity analysis, we tested the hypothesis that common underlying distributions for the Au-|-Au and 
Cu+Gu data can be used to represent the data. It is possible to also perform the same comparisons with the 
unidentified charged hadron spectra. In Fig. [TH a pair- fit using PHOBOS Cu-t-Gu data at .^snn — 200 GeV [s^l is 
invoked to extract the underlying mono- and multi-distributions, which are then used to reconstitute the -R^X'' 
Cu+Gu for each centrality bin. In contrast. Fig. 1151 uses the underlying mono- and multi-distributions derived from 
Au+Au collisions to form the i?^^ f'^'" each centrality bin. In both cases, a reasonable fit is found. For completeness, 
the right panel in Fig. [T3] shows a comparison between the underlying mono and multi distributions from Au+Au 
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FIG. 15. Nuclear modification factor (scaled by A'^part not A'coii) for unidentified charged hadrons from PHOBOS (filled band) 
in Cu+Cu collisions at y^SNN = 200 GeV using the model fit results for the mono and multi distributions from the Au+ Au 
analysis. The dashed (solid) lines use the results from fitting the PHOBOS (PHENIX) data. 



(symbols) and those derived from Cu+Cu (lines). The two sets of symbols for Cu+Cu represent the pair fit method 
using PHOBOS centrality bins 0-6% and 45-50% (open symbols) and 0-6% and 45-50% (closed symbols). 



FREEZE-OUT PROPERTIES 



As an extension to the low-pT particle studies, it has proven useful to fit the resultant spectra with models to derive 
freeze-out properties of the system. Commonly, the kinetic freeze-out temperature (Tkin) and radial flow velocity (/3) 
of the system are derived from a Blast-wave [35| fit to the low-p^ spectra of identified pions, kaons, and protons. For 



a description see for example Refs. [36|, |37|. At kinetic freeze-out all elastic collisions cease and the spectral shape 



is fixed. The extracted parameters reflect the state of the system at that time. What has been observed, through 
the Blast-wave fits, is that T^m and /3 vary quite strongly as a function of centrality. Here, we address the question 
whether those results can be described within the mono/multi framework. In a similar way to the above analysis, 
we consider that mono and multi have their own "universal" Tkin and /3 values. Using the least-x^ method described 
above, we fit the centrality dependence of the freeze-out variables and extract those parameters {Y variables in Eqn[T]) 
for mono and multi. The left and right panels of Fig. 1161 shows the result of fitting the freeze-out temperature (Tkin) 
and radial fiow velocity (/3), respectively, at -^snn = 200 and 62.4 GeV from STAR [s^. From the fits we find that 
the centrality dependence is well described in this framework, and the extracted Tkin and /3 values for the underlying 
mono and multi components are given in Table [Hi Although there is some collision energy dependence, the more 
striking feature is that the mono results for Tkin and /3 are significantly different from the multi values. A consistent 
picture emerges that the multi component has considerably larger ^ values and a freeze-out temperature of about 
one-half of that of the mono. 

Figure [iTl shows the results of a comparative analysis using the freeze-out parameters from Cu+Cu collisions at 
STAR [37[. Differences between Au+Au and Cu+Cu evident, especially at the highest energy. The heavier Au+Au 
system has a larger /? and smaller Tkin than those extracted from Cu+Cu collisions, see Fig. [T71 Although there are 
differences, we did not consider the systematic uncertainties in the data when creating the underlying distributions 
for the two collision systems, so it becomes difficult to draw any strong conclusions. In comparison to the freeze- 
out parameters extracted for -pp collisions, for Au+Au collisions at 200 GeV, the extracted Tki n fo r the underlying 
mono component is about 25% higher than that from the corresponding value (0.127±0.013_|36|). Similarly, the 
extracted radial fiow velocity is about 35% lower than the corresponding -p-p value (0.244+0.081 [36[). We should note 
here, that we do not expect the mono to be precisely as minimum bias pp. In particular, the extracted underlying 
mono distributions for multiplicity and charged hadron spectra are different from minimum bias pp, especially in the 
low-pT region. 

From this analysis, if we take the Blast-wave model at face value, it appears as though the kinetic freeze-out of 
the mono collisions, around the surface, occurs earlier than the multi, as implied by the different temperatures at 
freeze-out. This, in turn, suggests that the kinetic freeze-out occurs at a later time such that elastic collisions still 
occur in the multi after the complete freeze-out of the mono. The radial velocity, /?, is found to be far greater for the 
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multi, possibly indicating a more explosive expansion. Thus, one could argue that the multi system is simply more 
densely compact. We note that, the Blast-wave fits are simply reflecting the difference between the underlying mono 
and multi spectral shapes. 




FIG. 16. (color online) Least-x^ fit to the kinetic freeze-out temperature (left) and radial flow velocity (right) for Au+Au 
collisions at 200 (filled band) and 62.4 GeV (outline band). Dashed (dotted) lines represent the 200 GeV (62.4 GeV) fit results. 
The left (right) pointing arrows depict the mono (multi) Tkin or /3 in each panel. 




FIG. 17. (color online) Least-x'^ fit to the kinetic freeze-out temperature (left) and radial flow velocity (right) for Cu+Cu 
collisions at 200 (filled band) and 62.4 GeV (outline band). Dashed (dotted) lines represent the 200 GeV (62.4 GeV) fit results. 
The left (right) pointing arrows depict the mono (multi) Tkin//3. 

In the evolution of the collision, a chemical freeze-out occurs prior to the kinetic freeze-out. Chemical freeze-out 
occurs once all inelastic collisions have ceased; freezing the particle species. It is found that this temperature {Tch) 
is invariant with centrality [s^l, thus both the mono and multi must have the same apparent chemical freeze-out 
temperature. 

SUMMARY 

In summary, we have shown that some aspects of heavy ion data can be reproduced using a simple two-component 
model, with participant-like scaling variables. As an initial test of the model, we find that the multiplicity and 
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TABLE II. Mono and multi freeze-out parameters as derived from the least-x^ fit to the mid-rapidity data from STAR. 







Au-f Au 






Cu+Cu 




Energy 


mono 




multi 




mono 




multi 




(GeV) 


Tkin [GeV] 


P 


Tkin [GeV] 


P 


Tkin [GeV] 


/3 


Tkin [GeV] 




200 


0.160 


0.163 


0.084 


0.639 


0.152 


0.268 


0.112 


0.555 


62.4 


0.168 


0.184 


0.096 


0.576 


0.156 


0.226 


0.114 


0.548 



unidentified charged hadron spectra can be reproduced with an underlying gluon-dominated distribution; the latter 
without the need for a large suppression at high-pT- By fitting the data, we can simultaneously extract underlying 
distributions for each sub-component of the collision: the singly-hit mono participants and the multiply-hit multi. 
From this, a suppression is observed, however, it is limited to particles originating at surface-collisions (mono) and 
mostly in the low- to intermediate-pT region. The component characterizing the multiple nucleon-nucleon interactions 
is found to be peaked in the intermediate pt region and becomes flat at the highest measured pt- With these two 
sub-components, we observe that there is a mechanism by which one can reliably reproduce the multiplicity, spectral 
shapes, and freeze-out parameters. These is no need for a large (up to a factor of five) suppression which is currently 
used to describe the heavy ion data using the A'coii variable. At the center of this model picture is a dense (perhaps 
gluon dominated) system (multi) which is opaque to the particles produced at the periphery of the collision (mono) . 
A gluon-dominated picture is not excluded by other measurements, as discussed, the anomalous increase of baryons, 
with respect to mesons, could be expected from a system with a higher number of gluon-jets than quark-jets, though 
more details are outside the scope of the current paper. 

To conclude, in our model, we assert that the yield of high-px particles is not suppressed. The apparent reduction 
in yields observed at high-p-p is in fact due to a different particle production mechanism than what is expected from 
minimum bias pp interactions. In particular this mechanism does not produce high-p^ particles at the expected binary 
collision rate, but rather at a slower rate as determined by the geometry of the collision, in a similar way to that 
expected from close-packed gluons. 

This work was partially supported by US DOE Grants DE-FG02-04ER41325, DE-FG03-86ER40271, and DE-FG02- 
94ER40865. 
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